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1. Introduction 

The title of my talk is both ambitious and pretentious! I hasten to state that the mandate 
given to me is rather limited, namely to review the phenomenology of hadronic states discovered 
recently in the mass region of the charmonia and the bottomonia. Spearheaded by the experiments 
at the B factories and the Tevatron, with the experiments at the LHC as welcome new-comers, an 
impressive number of new states have been reported. Generically called X, Y and Z, these states 
defy a conventional quarkonia interpretation; this certainly holds for the majority of them. Their 
gross properties, such as the spin-parity assignments, masses, production mechanisms and decay 
modes, have been discussed in a number of comprehensive reviews [[I], g]. 

There have been a number of more recent developments in the field of quarkonium spec- 
troscopy and I will confine myself just to their discussion. They involve the observation of the two 
charged bottomonium-like resonances by the Belle Collaboration [||] in the 7T ± T(«5') (n = 1,2,3) 
and n ± hi,(mP) (m = 1,2) mass spectra that are produced in association with a single charged pion 
in e + e~ annihilation at energies near the Y(5S) resonance. Here h b {mP) are the P-wave spin-singlet 
bottomonia states. Calling the charged particles Z^(10610) and Z^(10650), their masses and the de- 
cay widths averaged over the five final states are, respectively, M[Zfo(10610)] = 10608.4±2.0 MeV, 
r[Z fc (10610)] = 15.6±2.5 MeV, and M[Z b ( 10650)] = 10653.2±1.5 MeV, r[Z fo (10650)] = 14.4 ± 
3.2 MeV. The favoured quantum number assignments for both are I G (J P ) = 1 + (1 + ). This discov- 
ery was preceded by the observation of the hb(lP) and h b (2P) states, also by the Belle Collabora- 
tion [@] in the reaction e + e~ -> h b (nP)K + %~ , with the masses M[h h (\P)} = (9898.25 ± 1.06+};^) 
MeV mdM[h b (lP)] = (10259.76 ±0.64+} ^) MeV. These measurements yield hyperfine splitting 
in the bottomonium sector, defined as the mass difference between the P-wave spin-singlet state 
hb(mP) and the weighted average of the corresponding P-wave triplet states, Xbj(nP), AMhf(«P) = 
(M(n 3 Pj)) -M(n l Pi), with AM HF (1P) = (1.62 ± 1.52) MeV and AM HF (2P) = (0.48+J;^) MeV. 
They are consistent with theoretical expectations and also with the hyperfine splitting measured in 
the charmonium sector AMhf = (0.14 ±0.30) MeV consistent with zero. Theoretically ex- 
pected widths of hb(lP) and h b (2P) are of order 100 keV [0], which are too small to be measured 
by Belle. 

Still on the subject of h b (lP), the BaBar collaboration [f7|] has presented evidence of its pro- 
duction in the decay Y(3S) — > 71%, (IP), followed by the decay h b (\P) — > yr]b(lS), in the distri- 
bution of the recoil mass against the 71° at the mass M[hb(lP)] = (9902 ±4± 1) MeV, which is 
consistent with the Belle measurements [Q]. The width of h},{\P) is consistent with the experimen- 
tal resolution, and the reported product branching ratio is 3§(Y(3S) — > n°hi,) x SS(h b — > yqb) = 
(3.7 ± 1.1 ±0.7) x 10~ 4 . In this, and also mM[h b (\P)}, the first error is statistical and the second 
systematic. The isospin- violating decay Y(3S) — > n Q h b {\P) is expected to have a branching fraction 
of about 10~ 3 [|, |], and the branching fraction m{h b {\P) -> yr] h {\S)) ~ (40-50)% [@1; hence, 
the measured product branching ratio is as anticipated theoretically. It is noteworthy that the decay 
Y(35) — > h b (\P)n + ii~ , which is suppressed by at least an order of magnitude compared to the 
decay Y(35) — > K Q h b (\P) [Sh, has not been observed. The observation of the singlet P-state in the 



charmonium sector h c (\P) has also been reported this year by the CLEO collaboration [10] in the 
process e + e~ — > 7i + 7i~h c: (lP) at the center-of-mass energy E CM _ = 4170 MeV. In fact, CLEO pio- 
neered the technique of searching for peaks in the mass spectrum recoiling against the 71°, and the 
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resulting mass M[h c (\P)} = (3525.27 ±0.17) MeV measured by this method is consistent with an 
earlier measurement of the h c {\P) mass from the decay \jf(2S) — > 7t°h c [11]. The product branch- 



ing ratio gg{\\f{2S) -> 7i°h c ) x &{h c -> yr\ c ) = (4.19 ±0.32 ±0.45) x 10~ 4 is in agreement with 
theoretical expectations, and is also very similar to what has been reported by Babar for the corre- 
sponding hb ( IP) product branching ratio, quoted above. However, there is an intriguing hint in the 
CLEO measurements of the cross section for e + e~ — > h c {\P)n + n~ , which rises at E CM . = 4260 
MeV. Since this is close to the mass of the J PC =\— hadron 7(4260), which is a candidate for the 
hidden cc tetraquark state, it would suggest that the mechanism e + e~ — > 7(4260) — > h c {\P)n + K~ 
has something to do with the rise in the cross section. This remains to be confirmed in the next 
round of precise experiments. 

2. Current experimental anomalies 

There is a number of anomalous features in the Belle data taken in the center-of-mass en- 



ergy region near the T(55) mass. The first of these was reported some three years ago [ [I4 |I3| ] 
in the processes e + e~ — > Y(lS)7t + 7i~ , Y(2S)n + 7i~ ,T(3<S)7T + 7r~, measured in the center-of-mass 
energy range between 10.83 GeV and 11.02 GeV. The enigmatic features of the Belle data are 

(i) the anomalously large decay widths (or cross sections) for the mentioned final states, and 

(ii) the dipion invariant mass distributions recoiling against the Y(1»S) and T(2S) states, which 
are at variance with similar spectra measured in the transitions involving lower mass bottomo- 
nium states T(nS) — > Y(mS)n + n~ (with m < n). To quantify the problem, the reported partial 
widths are r[r(15)w + w - )] = 0.59 ±0.04 ±0.09 MeV and T\I{2S)k + k-)\ = 0.85 ±0.07 ± 0.16 
MeV. Compared to the corresponding partial decay widths of the lower three T(nS) in = 2,3,4) 
states, r[T(25) -»■ Y(lS)n + n-)] ~ 6 keV, r[T(35) -> Y(2S)n + n-)} ~ 0.9 keV, and F[Y(4S) -»■ 
Y(lS)n + n~)] ~ 1.9 keV, the production of the Y(nS)7i + 7t~ in the energy region near the T(55) is 
larger by two to three orders of magnitude. The order keV partial widths are well-accounted for in 



the QCD multipole expansion [ |14] , |15[ ] based essentially on the Zweig-suppressed process shown 
in Fig. [I] (left-hand frame). The dipion invariant mass spectrum anticipated in the QCD multipole 
expansion is shown on the example of the decay T(4»S) — > Y(lS)n + 7i~ in Fig. [l] (right-hand frame) 
and compared with the data taken from the Belle collaboration at T(4»S) JT^]. They are in excellent 
agreement with each other. Not so, for the dipionic transitions measured in the T(55') region, in 
which the dipionic mass spectra are dominated by the scalar meson /o(980) and the tensor meson 
/ 2 (1270) (for the Y(lS)7i+7[- mode) and by the / (600) and / (980) mesons (for the Y(2S)x+7[- 
mode). This is illustrated in Fig. ^| for the process e + e~ — > Y(lS)n + n~ which shows the distribu- 
tions in the M n + n - (left-hand frame) and in the helicity angle (cos 6 distribution (right-hand frame). 
The dipion mass spectrum measured near the Y(5S) clearly shows peaks at /o(980) and /2(1270). 
An interpretation of the process in terms of the production and decay of a J PC = 1 tetraquark 
state [[T7|, [T8p (histograms and the solid curves) accounts well the experimental distributions. We 
will return to discuss the underlying dynamical model later in section 4 of this report. 

Not only are the cross sections for e + e~ — > Y(nS)7i + 7t~ (n = 1,2,3) near the T(55) anoma- 
lously large by at least two orders of magnitude, the same holds for the production of the P- 
wave spin-singlet bottomonia states h^imP) (m = 1,2), for which the production cross sections for 
e + e~ — > ht,(\P)n + n~ and e + e~ — > hi,(2P)n + n~ are also anomalously large [Hj. The ratios of the 
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Figure 1: Left frame: Zweig-suppressed diagram for the transition Y(nS) — > Y(mS)7zn with m < n, which 
forms the basis of the QCD estimates of the decay rates and distributions in heavy quarkonia dipionic tran- 
sitions. Right frame: The dipion invariant mass spectrum M K7l measured in the decay T(45) — > T(1S)JT7T by 
the Belle collaboration together with a theoretical curve based essentially on the diagram shown in the left 



frame. (From [ 16|.) 
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Figure 2: Fit results of the M n + n - distribution (a) and the cos0 distribution (b) for e + e~ — > Yj, — > 
Y(\S)n + n~ , normalized by the measured cross section by Belle [Q. The histograms represent theoret- 
ical fit results based on the tetraquarks hypothesis, while the crosses are the Belle data. The solid curves in 
(a) show purely resonant contributions from the /o(980) and /2(1270). (From 



production cross-sections in the indicated final states relative to that for the e + e 
production are as follows M] : 



T{2S)n+n- 



o[Y(\S)7l + K-} 

a[T(3S)7i + n-} 
a[h b {\P)% + %-\ 
a[h b (2P)iz + K^} 



0.638 ±0.065l°;°|? 

0.517 ±0.082 ±0.070 

0.407 ±0.07j£o7« 
0.78 ±0.091^ 



(2.1) 



We have already commented on the anomalous production cross sections in the Y(ns)Tt + TZ~ modes 
near the T(55) region. The ratios given in the last two equations above for the ht,(lP)n + n~ and 
hb(2P)7t + 7i~ are found to be of order unity, a feature which violates theoretical expectations as 
the processes T(55) — > ht,(mP)7i + 7i~ involve heavy quark spin-flip, which are suppressed by l/mj, 
in the amplitude. It is obvious that the production mechanisms of all five processes involving 
Y(nS)7i + ii~ (n = 1,2,3) and hb(mP)7t + 7i~ (m = 1,2) are exotic. In particular, the true mecha- 
nisms at work avoid the Zweig-suppression seen in similar dipionic transitions and evade power 
suppression due to the spin-flip transitions for the hb(mP)7i + 7t~ case. It is worth recalling that no 
excess of the kind seen in the Belle measurements near the T(55) Jl2| , [T3| , ^] is seen by them or 
any other experiment either at energies below or above the T(55) region. Any plausible theoretical 
explanation must account for all these features. 
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These measurements have invoked a number of theoretical ideas. Particularly interesting is the 



suggestion by Bondar et al. Q19[], in which the resonances Z^(10610) and Zfo(10650) are assumed 
mostly of a 'molecular' type due to their respective proximity with the B*B and B*B* thresholds. 
Thus, the internal dynamics of the states Z/,(10610) and Z/,(10650) is dominated by the coupling to 
meson pairs B*B — BB* and B*B*, respectively. In particular, the bb pair within the Z^(10610) and 
Zb( 10650) is an equal mixture of a spin-triplet and spin-singlet with the relative phase orthogonal 
between the two resonances, i.e., 

|Z,(10610)> = ^(<k®l&-lrt®°& 
ft (10650)) = ^ (0-^r aq + l; h ®0- Qq ) . (2.2) 



Here 0~ and 1 stand for the para- and ortho-states with negative parity. The assignments ( p.2| ) 
would predict that the mass difference M[Z fo (10650)] -M[Z fo (10610)] should be equal to that be- 
tween the B and B* masses. The observed mass difference of 46 MeV [ffl] is in neat agreement 



with this argument. The spin-structure in (22) also suggests that the resonances Z/,(10610) and 
Zb( 10650) have the same decay width. This again is in agreement within measurement errors with 
the Belle data [§: r[Zj,(10610)] = 15.6 ±2.5 MeV and T[Z b { 10650)] = 14.4 ± 3.2 MeV. The max- 
imal ortho-para mixing of the heavy quarks in the Zfo(10610) and Z/,(10650) resonances described 



by Eq. ( ]2.2[ ) also implies couplings of comparable strengths to channels with states of ortho- and 
para-bottomonium, leading to the following couplings of these resonances to the channels T(«5)7T ± 



and hbimP^im: 



C h E n T{nS) ■ (Zi(10610) -Z fc ( 10650)) , C T (p n x h b ) ■ (Z fe (10610) +Z fo (10650)) , (2.3) 

where Z/,(10610), Z/,(10650) and hb denote the polarization vectors of the corresponding spin-1 
states, and E K and p n are the pion energy and its three-momentum, respectively; C/, and Cr are a 



priori unknown coupling constants to be determined by data. The amplitudes described by Eq. ( p.3[ ) 
applied to the decays Y(5S) — > Y(nS)7i + 7t~ and Y(5S) — > hb(mS)n + 7i~ yield the right pattern of 
destructive and constructive interferences seen in the Dalitz distributions of these processes [Q]. All 
of these arguments are plausible. Further variations on the molecular theme and predictions can be 
seen in [21], g g]. 



However, the structure suggested in Eq. (2.2) is a postulate not yet seen in decays other than 



those of the Y(5S). A particular case in point are the decays of the Y(6S), where the available phase 
space for the decays Y(6S) — > Y(nS)n + n~ and Y(6S) — > hb(mP)n + n~ are much larger. Hence, the 
implications of Eqs. ( |2.2[ ) and ( |2.3[ ) should be, at least qualitatively, very similar to those discussed 
in the context of the Belle data from the Y(5S) region. This remains to be tested. In addition, 
there are also some specific features of the Belle data which do not go hand-in-hand with the usual 
understanding of a hadronic molecule, the closest example of which is the Deuteron. The masses 
of the Z/,(10610) and Z/,(10650) are above the respective thresholds. The Deuteron mass, on the 
other hand, lies below the threshold by about 2.2 MeV. Also, the decay widths of the Z^(10610) and 
Zb( 10650) are not particularly small, as one would expect for a hadron molecule. On the contrary, 
their decay widths are similar in order of magnitude as that of the Y(55). This is also curious as 
the other 'hadronic molecule' discussed at length in a similar context, namely the X(3872), has 
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a much smaller (by at least an order of magnitude) decay width, with the current 90% C.L. limit 
being T[X(3872)] < 1.2 MeV [g]. 

In the rest of this writeup, I will take the point of view that all the five anomalous processes 
measured by Belle at energies near the T(55) mass [12, [l3], have very little to do with the 
T(55) decays. Following Jl7| , [iq , [25] ], I will argue here that the final states Y(nS)n + n~ and 
hh(mP)7t + 7i~ are the decay products of the J PC = 1 tetraquark 7^(1 0890), which lies in mass 
tantalizingly close to the T(55) mass. More precise experiments are needed to tell the two apart 
than is the case currently. In the context of the Y(nS)ll + 7V~ final states, this was suggested in 
[17, 18, 25] and the dynamical model was shown to be consistent with the observed cross sections. 
Also, the measured dipion invariant mass distributions show the predicted scalar-and tensor-meson 
resonant structure. Moreover, in the tetraquark context, it is easier to understand why the production 
cross sections for e + e~ — > Yb( 10890) — > Y(nS)n + n~ , which involves a 3 P — > 3 5 transition, and for 
e + e~ — > 1^(10890) — > hb(mS)n + n~ , which involves a 3 P — > l P transition, are comparable to each 
other. Detailed distributions, including the resonant Z/j(10610) and Zfc(10650) effects are still being 
worked out in the tetraquark picture. 



3. Spectrum of bottom diquark-antidiquark states 

Much of the discussion of the tetraquark states involves the concept of diquarks (and anti- 
diquarks) as effective degrees of freedom, which will be used here to calculate the mass spec- 
tra, production and decay of the tetraquark states. In particular, four-quark configurations in 
the tetraquarks are assumed not to play a dominant role. Following this, the mass spectrum of 



tetraquarks [bq] [bq'] with q = u, d, s and c can be calculated using a Hamiltonian [ 26 ] 



H = 2m Q + H^ Q) +H^ Q) +H SL + H LL , (3.1) 



where: 



H SS Q) = 2 ( -Xbqh [fib ■ Sq) + (Sg ■ Sq)] , 

Hjs^ = 2(Jtbq)(Sb ' S«j + • S q ) + 2J^(Sfo • Sg) + 2J^g^(S q • Sq), 

H SL = 2A e (S i2 -L + S J -L), 

H LL = B Q LQQ{L f + l) . (3.2) 
All diquarks, denoted here by Q are assumed to be in the color triplet (3), as the diquarks in the 



(6) representation do not show binding [|27|]. Here thq is the constituent mass of the diquark [bq], 
{■3^bq)i is the spin-spin interaction between the quarks inside the diquarks, are the couplings 
ranging outside the diquark shells, Aq is the spin-orbit coupling of diquark and Bq corresponds 
to the contribution of the total angular momentum of the diquark-antidiquark system to its mass. 
The overall factor of 2 is used customarily in the literature. As the isospin-breaking effects are 
estimated to be of order 5-8 MeV for the tetraquarks p5| , [26] ], they are neglected in the 

mass estimates discussed below. 
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The parameters involved in the above Hamiltonian (3.2) can be obtained from the known 



meson and baryon masses by resorting to the constituent quark model Q29| ] 



ff = J>i + £2^(S r Sy), (3-3) 

i i<j 

where the sum runs over the hadron constituents. The coefficient Jfij depends on the flavour of the 
constituents i, j and on the particular colour state of the pair. The constituent quark masses and the 



couplings for the colour singlet and anti-triplet states are given in [|25|]. To calculate the spin- 
spin interaction of the QQ states explicitly, one uses the non-relativistic notation [28] \Sq,Sq; /), 
where Sq and Sq are the spin of diquark and antidiquark, respectively, and J is the total angular 
momentum. These states are then defined in terms of the direct product of the 2x2 matrices in 
spinor space, F a , which can be written in terms of the Pauli matrices as: 

!<> = ^;V = -Lo 2 o i , (3.4) 
V2 V2 



which then lead to the definition such as |0g,0g; 0/) = \ (02) ® (°2)- Others can be seen in [25]. 

The next step is the diagonalization of the Hamiltonian ( |3.1[ ) using the basis of states with 
definite diquark and antidiquark spin and total angular momentum., There are two different pos- 
sibilities [28]: Lowest lying states (Lqq = 0) and higher mass [bq][bq] states (Lqq = l). 
The [bq] [bq] states (Lqq = 0) can be classified in terms of the six possible states involving the good 
(spin-0) and bad (spin-1) diquarks (here, P is the parity and C the charge conjugation) 

i. Two states with J PC = 0++: 

|0 ++ > = |0 e ,0g;0 7 ); 

|0 ++, > = \IqAq; 0j). (3.5) 



ii. Three states with 7=1: 

|l ++ > = ^=(|0 fl ,l g ; l/) + |l fl) 0g; 1/)); 

|i +-/ > = i^>- (3.6) 

All these states have positive parity as both the good and bad diquarks have positive parity and 
Lqq = 0. The difference is in the charge conjugation quantum number, the state |1 ++ ) is even 
under charge conjugation, whereas 1 1 + ~) and 1 1 + ~') are odd. 

iii. One state with J PC = 2++: 

|2 ++ > = l g ; 2j). (3.7) 
Keeping in view that for Lqq = there is no spin-orbit and purely orbital term, the Hamiltonian 



(3.1) takes the form 



H = 2m Q + 2(^)3 [(S b ■ S q ) + (S s • S ? )] + 2J^ ? -(S ? • S 9 ) 

+2p^-)(S, • S q + S h • S q ) + 2je- bh (S b • S s ). (3.8) 
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The diagonalisation of the Hamiltonian (3.8) with the states defined above gives the eigenvalues 
which are needed to estimate the masses of these states. For the 1 ++ and 2 ++ states the Hamiltonian 



is diagonal with the eigenvalues [ ]28[ ] 



M (l ++ ) = 2m Q - {je bq )- 3 + X -X qq - - X bq + X -X bh (3.9) 
M (2++) = 2m [bq] + (X hq ) % + \jfT qq + X hq + \x bh . (3. 10) 



Mass of the constituent diquark can be estimated in one of two ways: We take the Belle data []12J as 
input and identify the 7^(10890) with the lightest of the 1 states, Fr^i, yielding a diquark mass 
m [bq] =5.251 GeV. This procedure is analogous to what was done in p8[], in which the mass of the 
diquark [cq] was fixed by using the mass of X (3872) as input, yielding m\ cq \ = 1.933 GeV. Instead, 
if we use this determination of my and use the formula mty = my + (m b — m c ), which has the 
virtue that the mass difference m c — m b is well determined, we get mi bq i = 5.267 GeV, yielding a 
difference of 16 MeV. This can be taken as an estimate of the theoretical error on m\u q \, which then 
yields an uncertainty of about 30 MeV in the estimates of the tetraquark masses from this source 
alone. For the corresponding ++ and l + ~ tetraquark states, there are two states each, and hence 
the Hamiltonian is not diagonal. After diagonalising the 2x2 matrices, the masses of these states 
are obtained. 

We now discuss orbital excitations with Lqq = 1 having both good and bad diquarks. Con- 
centrating on the 1 multiplet, we recall that there are eight tetraquark states (q = u,d), 
and the lightest isospin doublet is: 

M { '\ (S Q = 0,S Q = 0, S QQ = 0, L QQ = 1) = m [bq] + X, +B Q , (3.11) 



and the next in mass is: M$ (Sq = 1, Sq = 0, Sqq = 1, Lqq = l) = 2m^] + A + A2 — 2Ag +Bq, 
and so on. Values of = 1,2,3), Aq and fig are estimated in p5p. We identify the state 7^(10890) 



with My (in fact there are two of them, which differ in mass from each other by about 5-8 MeV, 
including isospin-breaking). This does not fix the quantity A, which is the mass difference of the 
good and the bad diquarks, i.e. A = mg (Sq = 1 ) — vuq (Sq = 0) . Following Jaffe and Wilczek [^7|] , 
the value of A for diquark [bq] is estimated as A = 202 MeV for q = u, d, s and c quarks. This is 
another source of potential uncertainty in estimating the tetraquark masses. The mass spectrum for 
the tetraquark states [bq] [bq] for q = u,d with J PC = ++ , 1 H h , 1 H , 1 and 2 ++ states is plotted 
in Fig. H in the isospin-symmetry limit. It is difficult to quote a theoretical error on the masses 
shown, with ±50 MeV presumably a good guess. Other estimates of the tetraquark mass spectra in 
the charm and bottom quark sectors can be seen in [|T|, 32, 33]. 



3.1 Estimates of the charged J p = 1 + tetraquark states 

In the tetraquark picture, one also anticipates a large number of charged states whose mass 
spectrum can be calculated in an analogous fashion as for their neutral counterparts just discussed. 
We would like to propose that the two charged J p = 1 + states Z^(10610) and Z^(10650) observed 
recently by the Belle Collaboration [S], and interpreted by them as the charged bottomonium states 
produced in the process T(5S) Z^(10610) + K T and T(55) -> Z^(10650) + 7i T , are indeed 
charged tetraquark states with the quark content = [bu] [bd] for the positively charged state (its 
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charge conjugate being Z b = [bu] [bd]). For the present discussion, they are produced in the decays 
of the J PC = 1 tetraquark 7^(1 0890). According to this interpretation, the decay chains involve 
F/,(10890) -> ((Z±(10610),Z± (10650) + jr* -> Y{ns)n + n-. A detailed dynamical model is un- 
der development with the aim of understanding the decay distributions in the kinematic variables 
available in these decays. 

We have estimated the masses of the isospin partners of Zfo(10610) and Zfo(10650), the two 
neutral 7 = 1 tetraquark states, denoted as 1 1 + ~) and 1 l + ~'). The 2x2 non-diagonal mass matrix 
for the neutral J PC = l + ~ states was, however, calculated numerically for A = 0. If we ignore the 
isospin-breaking effects in the tetraquark masses, which are small, then the charged counterparts 
have the masses M[Z & (10610)] = 10.386 GeV and M[Z 6 (10650)] = 10.527 GeV, given in Fig. |. 
As Zfe(10610) involves one good and one bad diquark and Z^(10650) involves two bad diquarks, 
including the A-dependent term, the non-diagonal 2x2 mass matrix gets modified to the following 
form 

M(l+-) = 2m Q + ~A- *<2±M + ( * ~ {Khqh + Kb " A K f'- Kh ~ b 
2 2 y K qq — K bh 2~^\ Kb ih~ 

The two eigenvalues can be written as E = ± x 2 + y 2 , with x = j + ( K"^)j — Ktq and y=Kqq — K bb , 
yielding 

M[Z,(10650)] = 2m Q + ~A - ^±^L + ^| + (^)- _ ^_)2 + {Kq _ _ ~f , (3 . 13) 

M[Z fe (10610)] = 2m Q + ^A - ^ ± *** - ^(| + (k^) 3 - K bq ) 2 + {x qq - K bb ) 2 . (3.14) 
Using the default values of the parameters | |25] ] 

m Q = 5.251 GeV, (k^) = 318 MeV, (K bl ) = 36 MeV, (k^) = 23 MeV, {K bq ) 3 =6MeV 

(3.15) 

we have now the following predictions for the two charged tetraquark masses 

M[Z fo (10610)] = 10.637 GeV; M[Z & (10650)] = 10.884 GeV, with A = 202 MeV . (3.16) 

These estimates are to be compared with the masses of the J p = 1+ states Z/,(10610) andZ,,(10650) 
reported by the Belle Collaboration [J3j] M[Z fo (10610)] = (10608 ±2.0) MeV and M[Z 6 (10650)] = 
( 10653.2 ± 1.5) MeV. They are in the right ball-park, but miss the measurements by approximately 
30 MeV and 230 MeV, respectively. More importantly, the mass difference between the two states 
has been measured precisely ^ M[Z h ( 10650)] -M[Z fo (10610)] ~ 45 MeV. The expression for this 
mass difference using the Hamiltonian ( ^^ ) is: 

M[Z fo (10650)] -M[Z b (10610)] = 2^J (| + (k^) 5 - K bq ) 2 + (K qq - K bb ) 2 . (3.17) 

The smallest value for the mass difference (140 MeV) is obtained for A = 0, which goes up to 247 
MeV for A = 202MeV. Both are larger than the measurements. Thus, the Belle data suggests that 
the Hamiltonian used here has to be augmented with an additional contribution. As the masses of 
the observed states Z 6 (10610) and Z 6 (10650) are rather close to the thresholds M(B) +M{B*) and 
2M(B*), respectively, this suggests that the threshold effects may impact on the masses and mass 
differences presented here. 
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10 528 fO"j " 10527(1+') 10 520(2 4 - + ) 

10504(1++) 1 ' I I 

10385 (0**) 10386(1+') 



0+ + 1+ + 1+ - 1 2 ++ 

Figure 3: Tetraquark mass spectrum with the valence quark content [bq] [bq] with q = u,d, assuming isospin 
symmetry. The value 10890 is an input for the lowest J c — I tetraquark state Fr^j. All masses are given 
inMeV. (From [||[.) 

4. Tetraquark-based analysis of the processes e + e — > Y(lS)(n + 7t~,K + K~,ri7t ) 

The cross sections and final state distributions for the processes e + e~ — > T( 15) (n + %~ , K + K , tj 7T°) 



near the Y(5S) have been presented in the tetraquark picture in [ |18| ] improving the results on the 
process e + e~ — > T(lS)7i + n published earlier The distributions for the process e + e~ — > 

T(2S)ti + 71 calculated in [ |r7j ] had a computational error, which has been corrected in the mean- 
while (see the Erratum in [|r7j]). These analyses are briefly reviewed in this section. Concentrating 
on the processes e + e~ — > Y (IS)(ti + 7i~ , K + K~ , T]7r°), there are essentially three important parts of 
the amplitude to be calculated consisting of the following: 

(i) Production mechanism of the J PC = 1 vector tetraquarks in e + e~ annihilation. To that 
end, we derive the equivalent of the Van-Royen-Weiskopf formula for the leptonic decay widths 
of the tetraquark states and Fr^i made up of a diquark and antidiquark, based on the diagram 
shown in Fig. |] (left-hand frame). 

24a 2 \Q[bu/bd]\ 2 



Yt 



r[\\o) 



(4.1) 



Here, Q\b u ] = 1/3 and Q\m = —2/3 are the electric charges of the constituent diquarks of the Y^ 
and Y\p^ , a is the fine-structure constant, the parameter K takes into account differing sizes of the 
tetraquarks compared to the standard bottomonia, with fc < 1 anticipated, and |/?j^(0)| 2 = 2.067 
GeV 5 [0] is the square of the derivative of the radial wave function for Xb(^P) taken at the origin. 
Hence, the leptonic widths of the tetraquark states are estimated as 



n?M -> e + e-) = 4r(7 N -> e+O « 83 K* eV , (4.2) 

which are substantially smaller than the leptonic width of the T(55) [^|. This is the reason why 
the states Fr M i and Y]p u ^ are not easily discernible in the /?/>-scan. Between the two, Yy bd -\ production 
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dominates and should be searched for in dedicated experiments. However, as the decays T(55) — > 
Y(nS)7i + 7i~ are Zweig-suppressed in the conventional Quarkonia descriptions, and hence have 
small branching ratios, the signal-to-background is much better for the discovery of the Yj,{ 10890) 
in the states Y {nS)n + tc~ . These, in fact, are the discovery channels of the 7^(10890) [£f|. 

(ii) The decay amplitudes for 7 fo (10890) -> Y(l^)(7U + n~ ,K + K~, 7771°) have non-resonant 
(continuum) contributions, as depicted in Fig. ^] (middle frame). They are parametrised in terms of 
two a priori unknown constants A and B , following JPfl ] : 



AC 



2 A B 3(9°) 2 jk?jk?-|q| 2 |k| 2 

•^o = T - r {h-k z )+ 1 2 1411 , 

JPJP' JPJP' is 

B\qm^ 

fpfp' s 

where the subscript denotes the 1 = part of the amplitudes, the superscripts 1C and 2C cor- 
respond to the 5- and D-wave continuum contributions, respectively, fp^ is the decay constant of 
pW, and |q|, k\ and are the magnitude of the three momentum of Yy and the energies of P and 
P' in the PP' rest frame, respectively. Using SU(3) symmetry results in the relations involving 
the various / = and / = 1 amplitudes: ^ C ' 2C (Y(\S)K+K-) = (V3/2) ^ C ' 2C (Y(IS)k + ti-), 
Jtl c - 2C {Y(\S)K+K-) =^ C ' 2C {Y(IS)K+K-) and ^^{Yi 15) tjw°) = s/2J?l c > zc {Y{\S)K+K-). 
We note that, in general, there is a third constant also present in the non-resonant amplitudes, char- 
acterising the term depending on the polarisation of the Yy. However, being suppressed by 1 /mi,, 
this is ignored. 

(iii) The resonant contributions, shown in the right-hand frame of Fig. f|, are expressed by the 
Breit-Wigner formula: 



yR _ 8RPP>gYlT(lS)R8e+e-Y b , (% 

Mp P , — m| + imgT R 



JK« = » V' ; h e m , (4.4) 



where I = for R = a, fo and f%, and / = 1 for 7? = Oq. The couplings for the scalar resonances 
S are defined through the Lagrangian Jz? = g SPP i(d ii P){d^P')S + gY b r(is)s^b^^S, while those for 
the f 2 are denned via «Sf = 2g hPpl (d ti P){d v P')f^ v + gY b r(is)f 2 Y b^vf2 V ■ The couplings g RPP < and 
8y't(\s)r nave mass dimensions —1 and 1, respectively. For the a, fo and a^, we adopt the Flatte 



model [35] and the details can be seen in 

With this input, a simultaneous fit to the binned Y(lS)7r + 7r~ data for the M n ^ n - and cos0 
distributions measured by Belle at y/s = 10.87 GeV JT2| ] were undertaken. Normalizing the dis- 
tributions by the measured cross section: dc n + n - jdM nn and dc n + n - j dcosQ , where o n + n - = 

CT r(is)^W CT r(1s)jr+jr- with a T( e i5)7t + ^ = 161 ±0 - 16 P b ©' tne results are shown in Fi §- I 
(histograms) and provide a good description of both the dipion mass spectrum and the angular 

distribution. 

The normalized M K + K - and M^o distributions are shown in Fig. || (a) and Fig. |5] (b), re- 
spectively. In these figures, the dotted (solid) curves show the dimeson invariant mass spectra 
from the resonant (total) contribution. Since these spectra are dominated by the scalars fo + c?q 
and a®, respectively, there is a strong correlation between the two cross sections. This is shown 
in Fig. H (c), where the normalized cross sections Ok+k- an d o^o are plotted resulting from 



the fits (dotted points) which all satisfy % /d.o.f. < 1.6 [18]. The current Belle measurement 
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Y, 





T(n5) e+ \ Yl , 
P 

e 

P 

Figure 4: Left frame: Van Royen-Weiskopf Diagram for the production of a J c = 1 tetraquark Y\ 
with the quark content [bu] [bu] in the process e + e~ —> y* — > Middle frame: Continuum contribution in 
the process e + e~ — > Y/, — > Y(nS)PP'. Right frame: Resonance contribution in the process e + e~ — > Yf, — > 
Y{nS)PP'. (Figures based on Q.) 





0.2 

&K+K- 

Figure 5: Predictions (a) of the M K+K - distribution for e + e~ -tY},-> Y(1S)K + K~ , (b) of the M^o distri- 
bution for e + e~ — s- Y(lS)r/7T° and (c) of the correlation between the cross sections of Y{\S)K + K~ and 
Y(lS)i~l7t°, normalized by the measured cross section for the Y(1S)k + k~ mode. In (a) and (b), the dotted 
(solid) curves show the dimeson invariant mass spectra from the resonant (total) contribution. In (c), the red 
dots represent predictions from the fit solutions satisfying ^ 2 /d.o.f. < 1.6. The shaded (green) band shows 
the current Belle measurement a K + K - = 0.1 1+°^ [01 . (From [jjj].) 



o K+K = 0.11j#o3 © is shown as a shaded (green) band on this figure. The tetraquark model [ |18| ] 



is in agreement with the Belle measurement, and prediction 1 .0 < CT^o < 2.0. will be further tested 
as and when the cross section a^o is measured. Another important test of the tetraquark model 



is [12] 



<*C(1S)K+K- = Q\bu\_ = 1 4 5 

Or(is)K°K° Q\ bd] 4' 

This remains to be tested. Finally, the corrected analysis [|l7|] of the dipion invariant mass spectrum 
and the helicity angle distribution (in cos0) for the process 7^(10890) — > Y(2S)7i + 7t~ are shown 
in Fig. ||, in which the normalization is given by the measured partial decay width T\Yb( 10890) — > 
Y(2S)ti + 71~} = 0.85±0.7±0.16 MeV [|P3|[. The dipion invariant mass spectrum is well accounted 
for also in this process (# 2 /d.o.f. = 12.6/7), but not the the angular distribution <ir/<icos 6. These 
distributions are being reevaluated taking into account the resonances Zfo(10610) and Zfo(10650). 

As a tentative summary of the tetraquark interpretation of the Belle data on e + e~ — > (Y(nS)n + 7t~ 
and e + e~ — > hb{mP)% + Tt~~ is that the existing analysis are encouraging and there exists a prima fa- 
cie case of its validity. However, the missing contributions from the charged tetraquarks in the 
analysis of the e + e~ — > (Y(nS)n + 7i~ data have to be incorporated and the fits of the e + e~ — > 
hb(mP)% + ii~ data have to be undertaken to get a definitive answer. 

I would like to thank Robert Fleischer and the organisers of the Beauty 201 1 conference for a 
very exciting meeting in Amsterdam. I also thank Christian Hambrock, Satoshi Mishima and Wei 
Wang for their help in preparing this talk and helpful discussions. 
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- [GeVl 




Figure 6: Dipion invariant mass (m K7l ) distribution (left-handed frame) and the cos0 distribution (right- 



handed frame) measured by the Belle collaboration for the final state T{2S)n + n |12| and the corresponding 
theoretical distributions (histograms) based on the tetraquark interpretation of the Yg( 10890). (From [|l7[|.) 
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